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Spectral imaging of the Sagittarius B2 region in multiple 7-mm 
molecular lines 
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ABSTRACT 

We have undertaken a spectral-line imaging survey of a 6 x 6 arcmin 2 area around Sgr B2 
near the centre of the Galaxy, in the range from 30 to 50 GHz, using the Mopra telescope. 
The spatial resolution varies from 1.0 to 1.4 arcmin and the spectral resolution from 1.6 to 
2.7 km s _1 over the frequency range. We present velocity-integrated emission images for 47 
lines: 38 molecular lines and 9 radio recombination lines. 

There are significant differences between the distributions of different molecules, in part 
due to spatial differences in chemical abundance across the complex. For example, HNCO 
and HOCO + are found preferentially in the north cloud, and CH2NH near Sgr B2 (N). Some 
of the differences between lines are due to excitation differences, as shown by the 36.17 and 
44.07 GHz lines of CH3OH, which have maser emission, compared to the 48.37 GHz line of 
CH3OH. Other major differences in integrated molecular line distribution are due to absorp- 
tion of the 7-mm free-free continuum emission (spatially traced by the radio recombination 
line emission) by cool intervening molecular material, causing a central dip in the molecular 
line distributions. 

These line distribution similarities and differences have been statistically described by 
principal component analysis (PCA), and interpreted in terms of simple Sgr B2 physical com- 
ponents of the cooler, lower density envelope, and dense, hot cores Sgr B2 (N), (M) and (S). 



Key words: ISM:individual (Sagittarius B2) 
ISM:kinematics and dynamics. 



ISM:molecules - radio lines:ISM 



I INTRODUCTION 

Sagittarius B2 (Sgr B2) is a giant molecular cloud complex near 
the centre of the Galaxy. There is spectacular star-formation ac- 
tivity, much of it deeply embedded in molecular cores, as traced 
by the high far-infrared lumin osity, compact and ultracompact H 

II regions (Gaumeetal. 1995), with maser emission from water 



tre spectral-line surveys dCummins etafl I 19861: iTurnerl 1 19891: 



iNummelin etafl 1998. 2000; Belloche et al. 2005, 200^). The PRI 



( McG rath. Goss & De Pree 2004), hy droxyl dGaume & Claussenl 
1 19901) . formaldehy de dMehringer et alj| 19941) and methanol, both 
class I and class II dCaswelllll996l : lMehringer & Mentenll 19971) . 

Three major star-forming centres are located in a north-south 
line, labelled north, middle and south, Sgr B2 (N), (M) and (S), 
which have strong continuum radio H II region free-free emission, 
and s ubmillimetre and millimetre wavelength th ermal dust emis- 
sion dGordon et a l. 1993; Pi erce-Price et alj2 000'). These cores are 
surrounded by a larger, less dense envelope. 

The hot cores Sgr B2 (N) and (M) are particularly rich 
in molecules, and so have been the targets of many millime- 
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MOS survey with the Greenbank Telescope (GBT) will survey 
Sgr B2 (N) over th e whole band from 300 MHz to 50 GHz 
dRemiian et all2008h . 

Sgr B2 is one of the most prominent features of the Cen- 
tral Molecular Zone (CMZ), the bar-shaped dSawada et ai1l2004l) 
molecular region in the central few hund red pc of the G alaxy, 
as shown in emission of ^ CO, 13 CO dOka et alj Il998l). C S 
dTsuboi. Handa & Ukitj[l999l) and HNCO dDahmen et alj|l997h . 
for example. 

The d i stance to Sgr B2 has been measured as 7.8lg .7 kpc by 
iReid et al.1 f2009t> . with the proper motion suggesting that Sgr B2 
is nearer by 0.13 kpc than the Galactic centre, and the projected 
distance from the Galactic centre is 0.09 kpc. For calculations, we 
assume a distance to Sgr B2 of 8.0 kpc. 

The Sgr B2 molecular cloud has complex kinematics, with the 
densest core centred on Sgr B2 (M) at velocity 60 - 65 km s , and 
this area correspon ding to a 'hole' in the emission in CO and CS 
at 40 - 50 km s" 1 dSato et alj|2000h . This has been interpreted by 
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lHasegawa et alj 11994. 2008]) as a cloud-cloud collision between a 
clump at 70 - 80 km s" 1 , and the 40 - 50 km s _1 cloud triggering 
the current star formation: the collision changes the velocity of ma- 
terial, leaving a deficit (the hole) in the 40 - 50 km s" 1 range, as 
material shifts to the intermediate velocity of the 60 - 65 km s^ 1 
range. 

The spatial distribution in the Sgr B2 molecular cloud is also 
complex, but can be largely described by a ridge of emission in 
a north-south line near the Sgr B2 (N), (M) and (S) h ot cores, or 
centred to the west of these cores Uones et al. 2008a). This ridge 
continues to a peak 1 arcmin north of Sgr B 2 (N), the 'north cloud' 
which is enhanced in HNCO and HOCO+ dMinh et al] 19981) . pos- 
sibly associated with the shock of the cloud-cloud collision. This 
ridge of emission coincides with the Sgr B2 extended env elope, as 
traced by sub-mm dust emission jPierce-Price et alj2000l) . 

In this paper, we present results of imaging the Sgr B2 region 
(6x6 arcmin 2 ) in multiple spectral lines, over 30 to 50 GHz (the 
7-mm band), using the Mopra telescope. This follows on a similar 
sp ectral-line imaging s urvey of Sgr B2 in the 3-mm band, prese nted 



jljones et al.1 d2008ah,|jones. Burton & Cunningham! d2008bh and 
Ijones. Burton & Lowej d2008ch . This study provides spatially re- 
solved line distributions of a large number of molecules and transi- 
tions to constrain chemical models of the Sgr B2 complex. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were made with the 22-m Mopra radio tele- 
scope, of the Australia Telescope National Facility (ATNF) using 
the UNSW-MOPS digital filterbank. The Mopra MMIC receiver 
has a bandwidth of 8 GHz, and the MOPS backend can cover the 
full 8-GHz range in the broad band mode. This gives four 2.2-GHz 
sub-bands each with 8192 channels of 0.27 MHz. There is also a 
zoom mode with up to sixteen spectra (up to four in each 2.2-GHz 
sub-band) of width 137 MHz. This provides higher spectral reso- 
lution with 4096 channels of 0.033 MHz in each zoom spectrum. 
We chose the broad band mode, as the lines in Sgr B2 are wide, 
so that the 0.27 MHz channels, corresponding to around 1.6 to 2.7 
km s , are quite adequate to resolve their structure. This mode 
also allowed a complete line-mapping survey to be made, without 
pre-selecting which lines would be covered, or the restrictions of 
selecting at most four lines in each 2.2 GHz sub-band. 

The Mopra receiver covers the range 30 to 50 GHz in the so- 
called 7-mm band, so the band actually covers wavelengths 6 to 10 
mm. We chose three tunings centred at 45.7, 38.0 and 34. 1 GHz, to 
cover the range whole 30 to 50 GHz range. The latter two tunings 
overlap, with two 2.2 GHz sub-bands in common in the range 34 to 
38 GHz. 

The area was observed with on-the-fly (OTF) mapping, in an 
sar Sg r B2 (N ) 
i et alj l2008ah . 

3-mm imaging of ljones et al. I fc008al) . we used galactic coordinates 
for the scan direction, and hence the square scan area. We used po- 
sition switching for bandpass calibration with the off-source refer- 
ence position observed before each 6 arcmin long source scan. The 
system temperature was calibrated with a noise diode. The 7-mm 
Mopra system does not use paddle scans (unlike the 3-mm system), 
and hence does not correct for the absorption through the atmo- 
sphere. This effect on calibration is around the 10 percent level, so 
we take that as the major uncertainty in the data. 

The OTF data were turned into FITS data cubes with the 



area 6x6 arcmin , cent red near Sg r B2 (N ) and (M), in a similar 
way to that described in Jon es et alj (2008a). However, unlike the 



LIVEDATA and GRIDZILLA packagefl The raw data files in RP- 
FITS^ format, were bandpass corrected using the off-source ref- 
erence spectra with LIVEDATA, a robust second order polynomial 
fitted to the base line and subtracted, and the data output as SDFITS 
(Garwood 2000) spectra. These spectra were then gridded into dat- 
acubes using GRIDZILLA, with a median filter for the interpolation. 
The median was used, as this is much more robust to the effect of 
bad data. 

The whole 8-GHz band for each tuning was first turned into 
four cubes with frequency pixels for the 2.2-GHz sub-bands, mak- 
ing twelve cubes overall covering the entire 30 to 50 GHz range. 
These were used to conduct a full line survey in this spectral range. 
We compiled a list of the strongest (identified) lines, as given in 
Tabled 

The GRIDZILLA scripts also allowed the lines (Table Q3 to be 
specfied, with their rest frequencies, so that the GRIDZILLA output 
provides FITS cubes with velocity coordinates. We used this for a 
second pass through the data, to produce individual cubes for each 
line, over a velocity range -300 to 300 km s" 1 . For the data over 
the frequency range 34 to 38 GHz, where two of the tunings over- 
lapped, we combined the data from both observations to improve 
the signal-to-noise. 

The FITS cubes were then read into the MIRIAD 
dSault. Teuben & Wright! 1 19951) package for further process- 
ing and analysis. As the emission is typically of low surface 
brightness, the data were smoothed in velocity, with a 3-point han- 
ning kernel, to make data cubes with improved surface brightness 
sensitivity. This gives around 0.54 MHz, or 3.2 to 5.4 km s" 1 
effective FWHM spectral resolution (across the spectral range), 
and ensures the data are Nyquist sampled in velocity. 

The resolution of the Mopra beam i n the 7-mm band was re- 
cently measured by Urquh art et al.l d201 Ol) to vary between 0.99 ar- 
cmin at 49 GHz and 1.37 arcmin at 31 GHz, or roughly 9 — X/D 
(where D is the 22-m dish diameter). We will assume these values 
to give resolution between 1.0 arcmin at 50 GHz, and 1.4 arcmin 
at 30 GHz, with the caveats that: a) this is somewhat smaller than 
previously quoted in the Mopra documentation as 82 arcsec (1.37 
arcmin) at 42 GHfi and; b) actually gives a somewhat flatter vari- 
ation than expected from 6 oc A. Since we are mostly concerned in 
this paper with the spatial and velocity structure, we have left the 
intensities in the T^ scale, without correction for the main beam 
efficiency onto the T m b scale. The measured main beam efficien- 
cies rjMB vary between 43 and 53 percent within the range 3 1 to 49 
GHz, and extended b eam efficiencies r\xs from 56 to 69 percent 
dUrquhart et~ai]|2010l) . 



3 RESULTS 

We have produced data cubes for the 47 strongest lines in the 30 
to 50 GHz range, as listed in Table Q] The flag in the last column 
of Table Q] refers to the source of the line identification. Most are 
from the NIST databas 43 of lines known in the interstellar medium 
dLovas &~Dr agoset 2004). As Sgr B2 is among the richest known 
sources of interstellar lines, these lines are mostly already well- 
known in Sgr B2. There are also some radio recombination lines, 



1 http://www.atnf.csiro.au/people/mcalabre/livedata.html 

2 http://www.atnf.csiro.au/computing/software/rpfits.html 

3 http://www.narrabri. atnf.csiro.au/mopra/obsinfo. html 

4 http://physics.nist.gov/PhysRefData/Micro/Html/contents.html 
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taken from the splatalogue comp ilatioiQ We have i dentified sev- 
eral lines in the JPL databastFl (Pickett et alj|l998h correspond- 
ing to other transitions of molecules (CH 2 NH, NH 2 CHO, NH 2 CN, 
CH3CHO) kn own in Sgr B2 from other frequencies, notably in the 
3-mm band in ljones et alj J2008ah . 

The RMS noise in the (Hanning smoothed) data cubes varies 
from 34 to 76 mK (T^). 

There are also a few other weaker lines detected (mostly at the 
Sgr B2 (N) and (M) cores), but as they are quite weak and we do 
not have confident line identifications, we not consider them further 
here. Other projects are in progress to obtain line survey data for 
Sgr B2, notably the PRIMOS^ survey with the 100-m Green Bank 
Telescope (GBT) of Sgr B2 (N) which will cover this 7-mm band 
dRemiian et ai] |2008l) . These line surveys have longer integration 
time and higher sensitivity, and are better for studying the weak 
lines in the cores than the mapping observations presented here. 

To study the spatial distribution of the different lines, we have 
made integrated emission images (Figs. Q] and [3] to [6j, We inte- 
grated each data cube over a velocity range which included signifi- 
cant emission (above around 3 a in the cubes) for each line, which 
means the velocity range is different for each line. 

The integrated images in Figs.[T]and[3]to[6]are plotted as grey- 
scale, with the scale bar to the right in K km s - 1 (TJ). The contour 
levels are in equal linear steps, mostly 2 K km s _1 , but sometimes 
5 or 10 K km s" 1 , and 100 K km s _1 for the very strong CH 3 OH 
maser line at 36.17 GHz. We plot as fiducial marks the positions 
of radio sources as crosses and mid-IR sources as squares, as in 
Ijones et al] d2008al) . As the axes of Figs. Q] to [6] are Galactic coor- 
dinates, the north-south line of cores Sgr B2 (N), (M) and (S) is 
at a (Galactic) position angle around 45 degrees, with Sgr B2 (N) 
being the cross (radio source) near the centre at I — 0.680 deg, 
b — —0.028 deg, Sgr B2 (M) being the cross and square (radio 
and infrared source) below and to the right at I = 0.667 deg, 
b = —0.036 deg and Sgr B2 (S) the cross and square further below 
and to the right at I = 0.657 deg, b = —0.041 deg. The radio and 
infrared peaks are plotted (with the N, M and S labels) in Fig. lAll 
in equatorial coordinates. 

We discuss these images in more detail below, in subsections 
l3~T1to l3~T4l 

For quantitative analysis of the differences in spatial distribu- 
tion of the different lines, we have fitted the position of the peak 
emission in the data cubes, with the MIRIAD task maxfit, which fits 
the position with a parabolic fit to the spatial pixels arou nd the peak 
pixel in the cube. This is similar to the analysis used in I Jones et al.l 
(2008a) for the 3-mm Sgr B2 Mopra imaging, but we only fit the 
strongest peak, not multiple peaks, as the resolution is lower here at 
7 mm and much of the spatial structure is merged together. As the 
emission is typically complex and extended, we do not generally fit 
the spatial structure as gaussian peaks. 

These peak positions are tabulated in Table [2] arranged in 
groups of similar lines, as plotted in Figs. Q] and [3] to [6] We note 
that due to absorption in some of the spectra, and line width differ- 
ences across the area, that the position of the peak emission in the 
data cubes is not the same as the position of the peak of integrated 
emission. 

The velocity structure was studied by fitting the spectrum in 
each data cube at the position of the pixel of strongest emission, 



Table 1. Summary of strong lines detected in the Mopra observations. The 
first column gives an approximate frequency we have used for labelling 
transitions in this paper. The second column identifies the species, and the 
third its transition. The last columns indicates the rest frequency and the 
database used for the identification: L, Lovas: S, splatalogue; J, JPL (see 
text). The rest frequencies with a * indicate the frequency used for lines 
corresponding to multiple transitions, most also indicated as group (gp) in 
the transition list. 



Approx. 


Line ID 




Exact 




Freq. 


molecule 


Transition 


Rest Freq. 




GHz 


or atom 




MHz 




30.00 


SO 


not - 0( 11 


30001.547 


L 


31.22 


H 


59 a 


31223.313 


s 


31.95 




12-11 


31951 777 

j 1 y j l . / / / 


L 


32.85 


[r 


58 a 


32852 196 

J^OJfa. L y\J 


s 


33.75 


CCS 


3,2 - 2,1 


33751.370 


L 


34.60 


[r 


57 a 


34596 383 


s 


34.61 




13-12 


14614 185 


L 


35.07 


CH 2 NH 


1(0 31 - 2(1 21 pn 


35065 545* 


j 


35.27 


H 13 CCCN 


4 — 3 vnriiin 


35267.440* 


L 


36.17 


CH3OH 


4(— 1,4) - 3(0,3) E 


36169 290 

JUt .i. yu 


L 


36 39 




*+ — J glUlip 




L 


36 47 




JO UL 


jUtOO.ZO 




36 49 


OCS 


3 — 2 


JUtOO.O 1 J 


L 


36 80 


CH3CN 


9ffU — 1 f(Tl an 




L 


37 28 


i-trvN 


14—13 


j 1 z. 1 \j. yo j 


L 


38 47 


T-J 


55 a 


JO^ / J . J JO 


3 


38 506 


V. 1 V I Ivy 


2(0 2 1 - 1 (0 n F 


38506 035 


L 


38.512 


fH., PHO 


2(0 21 - 1(0 11 A 


385 1 2 08 1 


L 


1Q 16 


Pt-T., PHPi 
v n \ \ [iv' 


9 n 1 ^ 1 n n\ p 

Z.\ 1,1) — LyL^KJl d 


J7JOZ.JJJ 


T 


J7.J7 


PT-T., PT-T Pi 
v n . ; v [ivy 


9C1 1 ^ 1 n Oi\ A 




T 


39 73 


Nt-T^PN 




J7 / Z.J.JOI 1 


J 


1Q Q4 




15 — 14 


jyyDy.J /"+ 


L 


40 25 


CH 2 CN 


J 1 rri'oi 1 W 
_ — 1 glULip 


4.09 884* 

tUZ IJ.OOi 


L 


40.63 




54 a 


HUUJU.tyO 


3 


40.88 


I'lrT^v v i [ v v 


9C1 7\ 1C1 1~1 on 


4.0S7S llf\fr* 
tuo / j.z. / 


j 


42 39 


Nt-ToPt-TO 
[vn2Lnu 


9CO 9^ Iff) n<m 


tZJOO.U / u 


j 


42 60 




16—15 


49(^09 1 


L 


42 67 


HCS^ 


1—0 


AlfHA 1 Q7 


L 


42 77 


HOCO+ 


9C0 9*1 1 CO 1^1 


tz / . 1 y 1 j 


L 


42.88 


29 SiO 


1-0 v=0 


42879 922 


L 


42.95 


[r 


53 a 


42951 968 

*T^y *J L .7UO 


s 


43 42 


SiO 


1—0 v=0 


A^Al^, SAzt 
t-Jt-Zj . OD4 


L 


43.96 


HNCO 


2(0,2) - 1(0,1) gp 


43962.998* 


L 


44.07 


CH3OH 


7(0,7) - 6(1,6) A+ 


44069.476 


L 


44.08 


H 13 CCCN 


5-4 


44084.172 


L 


45.26 


HC 5 N 


17-16 


45264.720 


L 


45.30 


HC 13 CCN 


5-4 


45297.346 


L 




HCC 13 CN 


5-4 


45301.707* 


L 


45.38 


CCS 


4,3-3,2 


45379.029 


L 


45.45 


H 


52 a 


45453.719 


s 


45.49 


HC 3 N 


5-4 group 


45490.316* 


L 


46.25 


13 cs 


1-0 


46247.580 


L 


47.93 


HC5N 


18-17 


47927.275 


L 


48.15 


H 


51 a 


48153.597 


s 


48.21 


C 34 S 


1-0 


48206.946 


L 


48.37 


CH3OH 


1(0,1) -0(0,0) A+ 


48372.467* 


L 




CH3OH 


1(0,1)- 0(0,0) E 


48376.889 


L 


48.65 


OCS 


4-3 


48651.6043 


L 


48.99 


cs 


1-0 


48990.957 


L 



5 http://www.splatalogue.net/ 

B http://spec.jpl.nasa.gov/ftp/pub/catalog/catform.html 
7 http://www.cv.nrao.edu/ aremijan/PRIMOS 
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Table 2. Summary of peaks fitted to the lines. The lines are grouped into the 
same order as plotted in Figs.[T]and[3]to[6] The positions are the fits to the 
maximum in the channel of strongest emission, and the velocity centre, ve- 
locity width (FWHM) and peak are from gaussian fits at the brightest pixel 
of the line emission. Note that the velocity centres use the rest frequencies 
from Table [T] so will be confused for the groups of lines, and the velocity 
widths will be increased for such confused lines. 



Line ID 
molecule 
or atom 



Rest 
Freq. 
GHz 



Lat. 

peak 
deg 



Long, 
peak 
deg 



Veloc. Veloc. 



Peak 



centre width T"5 
kms" 1 K 



H59« 


31.22 





675 


-0 


031 


62 


48 


0.40 


Ft 58 a 


32.85 





674 


-0 


034 


65 


38 


0.45 


Ft 57 a 


34.60 





672 


-0 


035 


64 


34 


0.46 


Ft 56 a 


36.47 





676 


-0 


034 


65 


36 


0.40 


Ft 55 a 


38.47 





674 


-0 


035 


65 


34 


0.43 


H54a 


40.63 





672 


-0 


036 


62 


46 


0.40 


H53 a 


42.95 





670 


-0 


034 


67 


35 


0.33 


Ft 52 a 


45.45 





672 


-0 


034 


63 


32 


0.33 


H51 a 


48.15 





675 


-0 


034 


64 


54 


0.39 


HC 3 N 


36.39 





692 


-0 


021 


65 


31 


2.58 


HC 3 N 


45.49 





691 


-0 


022 


66 


30 


2.17 


HCC 13 CN a 


45.30 





682 


-0 


018 


84 


43 


0.17 


H CCCJN 


35.27 





694 


-0 


010 


65 


25 


0.20 


H 13 CCCN 


44.08 





679 


-0 


019 


69 


26 


0.19 


HC 5 N 


31.95 





699 


-0 


022 


65 


24 


0.29 


HC 5 N 


34.61 





691 


-0 


015 


65 


21 


0.33 


HC 5 N 


37.28 





685 


-0 


010 


67 


22 


0.31 


HC 5 N 


39.94 





675 


-0 


015 


68 


21 


0.34 


HC 5 N 


42.60 





688 


-0 


021 


69 


20 


0.32 


HC 5 N 


45.26 





698 


-0 


024 


64 


21 


0.27 


HC 5 N 


47.93 





696 


-0 


025 


65 


25 


0.39 


CH 3 OH 


36.17 





693 


-0 


024 


67 


22 


35.7 


CH3OH 


44.07 





672 


-0 


028 


65 


21 


1.77 


CH3OH 


48.37 





689 


-0 


022 


64 


35 


2.59 


CH3CHO 


38.506 





687 


-0 


014 


70 


26 


0.26 


CH3CHO 


38.512 





679 


-0 


007 


64 


39 


0.22 


CH3CHO 


39.36 





687 


-0 


021 


70 


18 


0.16 


CH3CHO 


39.59 





684 


-0 


012 


68 


20 


0.21 


NH 2 CHO 


40.88 





684 


-0 


008 


64 


18 


0.21 


NH2CHO 


42.39 





681 


-0 


012 


70 


25 


0.24 


SO b 


30.00 





695 


-0 


008 


58 


(25) 


0.26 


SiO 


43.42 





707 


-0 


021 


57 


50 


0.67 


29 SiO 


42.88 





673 


-0 


009 


67 


25 


0.14 


CS c 


48.99 





664 


-0 


040 


51 


11 


1.91 














83 


26 


1.39 


13 cs 


46.25 





656 


-0 


041 


60 


25 


0.31 


c 34 s 


48.21 





662 


-0 


042 


61 


33 


0.45 


CCS 


33.75 





693 


-0 


020 


66 


24 


0.27 


CCS 


45.38 





683 


-0 


008 


65 


18 


0.24 


CH3CN 


36.80 





680 


-0 


012 


72 


33 


0.71 


ocs 


36.49 





676 


-0 


032 


63 


20 


0.47 


ocs 


48.65 





692 


-0 


024 


62 


25 


0.52 


HNCO 


43.96 





696 


-0 


021 


64 


24 


2.30 


HOCO+ 


42.77 





682 


-0 


011 


67 


20 


0.34 


CH 2 NH 


35.07 





668 


-0 


034 


60 


20 


0.48 


NH 2 CN 


39.73 





699 


-0 


024 


67 


47 


0.20 


HCS+ 


42.67 





674 


-0 


011 


66 


19 


0.17 


CH 2 CN 


40.25 





651 


-0 


055 


54 


69 


0.21 



a line labelled as HCC 13 CN a blend with HC 13 CCN so larger line width 
and biased velocity 

b SO line at band edge, fitted with fixed width, due to poor data 
c CS line with self-absorption, fitted as two gaussians 




0.69 0.68 0.67 0.66 

GL0N (degrees) 

Figure 2. The fitted peak positions of the 7-mm radio recombination lines 
(filled) triangles, relative to the high resolution 7-mm radio continuum emis- 
sion, from ATCA observations, in grey-scale and contours, with the low 
resolution Mopra integrated Ft 51 a also plotted as contours. The recombi- 
nation lines include emission from Sgr B2 (N), centre, Sgr B2 (M), lower 
right, and the extended envelope which is resolved out in the ATCA data. 



with the MIRIAD task gaufit. The fitted velocity and velocity width 
are given in Table [2] The peak brightness (TX) at this position is 
also listed, to give an indication of line strength. We have used the 
hanning smoothed data cubes, with spectral resolution 3.2 to 5.4 
km s~ 1 , for the fitting, but this is not expected to significantly affect 
the velocity widths. 



3.1 Hydrogen recombination lines 

There are nine radio recombination lines (RRLs) of hydrogen, in 
the series H 51 a to H 59 a, shown in Fig. [7] They peak between 
the two strong free-free radio continuum sources Sgr B2 (N) and 
(M), with both sources contributing to the RRL emission. The mean 
peak position for the nine lines (Table O is I — 0.673 deg, b = 
—0.034 deg, mean velocity 64 km s" 1 and mean velocity width 
40kms~\ 

We show in Fig. [2] the peak positions of the RRLs relative 
to (our previously unpublished) higher resolution Australia Tele- 
scope Compact Array (ATCA) observations of the continuum. The 
low resolution Mopra data merge together the RRL emission from 
Sgr B2 (N) and (M). We also note (see the Appendix) that the 7-mm 
free-free emission has an extended diffuse component, in addition 
to the cores Sgr B2 (N) and (M), which is resolved out in the ATCA 
observations. The extended component contributes close to half of 
the radio flux in the 30 to 50 GHz range. 

For comparison BIMA observ ations of H 59 a RRL in both 
Sgr B2 (M) and (N) are given in |Pei. Liu & SnvdeJ d2000h and 
higher resol ution VLA observati ons of the H 52 a in Sgr B2 (M) 
are given bv lde Pree et al. (1996). 

The 7-mm free-free emission probed by the RRLs and shown 
in Fig.[2]also indicates where there may be absorption of molecular 
lines (Figs.[3]to[6]l by the strong Sgr B2 (N) and (M) continuum. 
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Figure 1. The integrated emission images for the nine recombination lines of hydrogen H 59 a to H 51 a. In this figure and later images the crosses (+) denote 
the position of radio peaks and the open squares (□) positions of infrared peaks. The line of cores Sgr B2 (N), (M) and (S) is at a (Galactic) position angle 
around 45 degrees, with Sgr B2 (N) being the cross near the centre at I = 0.680 deg, b = —0.028 deg, Sgr B2 (M) being the cross and square below and to 
the right at I = 0.667 deg, b = —0.036 deg and Sgr B2 (S) the cross and square further below and to the right at I = 0.657 deg, b = —0.041 deg. For clarity, 
these cores are labelled in this figure, and later figures. The scale bar is in units of K km s_i, for J Tjjdu. 

3.2 HC 3 N and 13 C isotopologues 

There are two lines of cyanoacetylene HC3N, shown in Fig. [3] the 
4-3 (36.39 GHz) and 5-4 (45.49 GHz) groups of transitions. The 
integrated emission is distributed on a ridge at higher latitudes than 
the Sgr B2 (N), (M) and (S) line (shown as crosses), that is to the 
north-west in equatorial coordinates, with the peak to the north of 
Sgr B2 (N). This is similar to the distributi on of the HC3N 3-m m 
lines 9-8, 10-9, 11-10 and 12-11 found in I Jones et al] J2008ah. It 
is also simil ar to the distribution of the 4- 3 line inllis & Goldsmith! 
dl99ll) and IChung. Ohishi & Morirriotcl i 19941) . but these latter 
data covered smaller areas. 

There are also three lines of the 13 C isotopologues of HC3N, 
the combined HCC 13 CN and HCC 13 CN 5-4 line (45.30 GHz), the 
H 13 CCCN 4-3 (35.27 GHz) and H 13 CCCN 5-4 (44.08 GHz) lines. 

© 2010 RAS, MNRAS 0Q0.[THT8l 



These are all quite weak lines, but the 45.30 GHz line shows good 
agreement in spatial distribution to the HC3N lines. The other two 
show fair agreement, but may be affected by scanning artifacts, 
with stripes along the longitude direction. 

The mean peak position for the five lines is I = 0.688 deg, 
b — —0.018 deg, and, for the four lines not blended, mean velocity 
66 km s _1 and mean velocity width 28 km s _1 . (Note that the 
blended 45.30 GHz line has larger width, and anomalous velocity, 
in TableU) 



3.3 HC 5 N 

There are seven lines of cyanodiacetylene HC5N, shown in Fig.[3] 
the series from 12-11 (31.95 GHz) to 18-17 (47.93 GHz) transi- 



P. A. Jones et al. 




Figure 3. The integrated emission images for two lines of HC3N, three lines of 13 C isotopologues of HC3N, and seven lines of HC5N. 
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tions. The lines are all quite weak, but the distribution is similar to 
that of HC 3 N (above), best seen in the 31.95 GHz and 39.94 GHz 
lines, but broadly similar in the other lines, with the ridge to more 
positive latitudes than the radio sources (crosses). The mean peak 
position for the seven lines is I = 0.690 deg, b — —0.019 deg, 
mean velocity 66 km s _1 and mean velocity width 22 km s -1 . 



3.4 CH 3 OH 

There are three lines of methanol CH3OH, shown in Fig. [4] the 
4(-l,4) - 3(0,3) E (36.17 GHz), 7(0,7) - 6(1,6) A+ (44.07 GHz) 
and the 1(0,1) - 0(0,0) A+ plus 1(0,1) -0(0,0) E (48.37 GHz) lines. 
The three line have quite different distributions, due to the complex 
collisional and radiative excitation processes of methanol. 

The 36.17 GHz line can act as a (class I) maser, which gives 
the high brightness p eak seen here in F ig. [4] Note that the peak 
is at the north cloud J Jones et aT]|2008 j) at I = 0.693 deg, b = 
— 0.024 deg. The dis t ributio n is in good agreement with that found 
bv lLiechti & Wilsor] d 19961) for this line, with higher resolution Ef- 
felsberg 100-m telescope observations. 

The 44.07 GHz line can also act as a (class I) maser. 
lMehringer & Menten j 19971) have used the VLA at much higher 
resolution than the Mopra data here to map the distribution of 
this line. They find many maser and quasi-thermal peaks near the 
Sgr B2 (N) and (M) cores. The Mopra data (Fig.f4) are consistent 
with this distribution. 

The 48.37 GHz line has a distribution that is simil ar to HC3N 
(Fig. [3 ]> and other 3-mm transitions of methanol in Ijones et all 
d2008al) . This is not a m aser transition so co mparison with the 44.07 
GHz line is instructive dVal'ttsetal.lll99ll) . 



3.5 CH3CHO 

There are four lines of acetaldehyde CH3CHO, shown in Fig. [4] the 
2(0,2) - 1(0,1) E (38.506 GHz), 2(0,2) - 1(0,1) A (38.512 GHz), 
2(1,1) - 1(1,0) E (39.36 GHz) and 2(1,1) - 1(1,0) A (39.59 GHz) 
lines. 

The distribution of integr ated emission from aceta ldehyde 
is widespread, as found by IChengalur & Kanekarl ( 120031) . How- 
ever, we find for these 7-mm lines that there is a deficit of 
emission near Sgr B2 (N) a nd (M), probably due to absorption. 
IChengalur & KanekaJ J2003h find the 1(1,1) - 1(1,0) transition at 
1065 MHz in emission near Sgr B2 (N), (M) and (S), which they 
attribute to a weak maser. The distribution of the four lines is simi- 
lar, although the 39.36 and 39.59 GHz lines are weak. The emission 
peaks at the north cloud with mean peak position I = 0.684 deg, 
b = —0.014 deg, mean velocity 68 km s _1 and mean velocity 
width 26 km s , 



3.6 NH2CHO 

There are two lines of formamide NH2CHO, shown in Fig. [4] the 
2(1,2) - 1(1,1) group (40.88 GHz) and the 2(0,2) - 1(0,1) group 
(42.39 GHz). The distribution of integrated emission is similar 
to that of acetaldehyde (above) with widespread emission, peak- 
ing at the north cloud (mean peak position I — 0.683 deg, b — 
-0.010 deg) and a deficit near Sgr B2 (N) and (M). This distribu- 
tion is somewhat diff erent to that of the formamide 5(1,4) - 4(1,4) 
line at 102.07 GHz in I Jones etalj|2008al) which peaks on the ridge 



to the west of Sgr B2 (M), but this line is confused with the pro- 
tonate d formaldehyde H 2 COH+ 4(0,4) - 3(1,3) line dOhishi et al.l 
1 1996b . 



3.7 SO 

There is one line of sulphur monoxide SO, shown in Fig. [5] the 
1(0) - 0(1) (30.00 GHz) line. There is widespread emission over 
the area, with strong absorption near Sgr B2 (M). This is in con- 
trast to the 3-mm SO lines, where the emission is strongly peaked 
at Sgr B2 (M) in the 2(2) - 1(1) (86.09 GHz) and 2(3) - 1(2) 
(109.25 GHz) lines, and extended along the ridg e-line to the west of 
Sgr B2 (M) in the 3(2) - 2(1) (99.30 GHz) line dJones et alj2008at 
iGoldsmithetalJ 19871) . 



3.8 SiOand 29 SiO 

There are two lines of silicon monoxide SiO 1-0 v=0 (43.42 
GHz), and the isotopologue 29 SiO 1 - v=0 (42.88 GHz) shown in 
Fig. [5] The integrated emission is widely distributed in both lines, 
with absorption near Sgr B2 (M) and (N), like the SO line (above). 
The ratio between the isotopologue 29 SiO and SiO, indicates that 
the SiO is optically thick. This leads to the differences in the in- 
tegrated distribution, and fitted mean velocity and velocity width 
in Table [2] as the suppression of emission in the SiO line centre 
shifts the mean and increases the fitted width. The SiO 1-0 line 
distrib ution coincides with the pea k labelled SiO+69-0.06 of this 
line in iMartin-Pintado et al.l dl997l) . and has similar absorption at 
Sgr B2 (M) and (N) seen in the 3-mm SiO 2 - 1 (86.85 GHz) line 
in that p aper. I t is als o si milar to the distribu tion of the SiO 2 - 
1 line in iMinhl d2007l) and I Jones et all {2008a), but the latter is at 
lower signal to noise. 

We also have data for the maser lines SiO 1 - v=l (43.122 
GHz) and 1 - v=2 (42.821 GHz) but did not detect these lines 
above a 3 a level of 0. 10 K, although these line s have been detected 
in Sgr B2 by, for example Zapata et al. (20 0J^). If the v=l line was 
a strong as that reported by Zapata et al. ( 2009), corresponding to 
0. 18 K, it would have been detected here, but such masers are vari- 
able, so it is quite plausible that the strength has changed. 



3.9 CS, 13 CS and C 34 S 

There are three lines of carbon monosulphide CS 1 - (48.99 GHz) 
and the isotopologues 13 CS 1 - (46.25 GHz) and C 34 S 1 - 
(48.21 GHz) shown in Fig. [5] Like the SiO line above, there is 
widespread emission, with absorption near Sgr B2 (M) and (N). 
There are optical depth effects which make the integrated emission 
of the CS line and the isotopologues 13 CS 1 - and C 34 S differ- 
ent. The peaks in the data cubes (Table |2]l are near Sgr B2 (M) at 
peak position I — 0.661 deg, b = —0.041 deg (mean of the five 
lines). The CS line shows absorption, so is fitted in Table|2]as two 
gaussians, offset in velocity. The peaks of integrated emission (Fig. 
[5} are at different positions, due to the effe cts of absorption. T hese 
distributions are similar to that found by I Jones et alT d 2008al) for 
the 3-mm lines CS 2 - 1 (97.98 GHz) and the isotopologues 13 CS 
2 - 1 (92.49 GHz) and C 34 S 2 - 1 (96.4 1 GHz). The CS 1 - 
(48.99 GHz) line has also been mapped bv lSato et al] d2000h . with 
the Nobeyama 45-m telescope at higher resolution, showing similar 
line distribution. 
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Figure 4. The integrated emission images for three lines of CH3OH (methanol), four lines of CH3CHO and two lines of NH2CHO. 



3.10 CCS 

There are two lines of dicarbon monosulphide CCS shown in Fig. 
|5] the 3,2 - 2,1 (33.75 GHz) and 4,3 - 3,2 (45.38 GHz) lines. 
The emission peaks in the north cloud, with mean peak position 
/ = 0.688 deg, b = —0.014 deg, mean velocity 66 km s _1 and 
mean velocity width 21 km s _1 . The integrated emission avoids 
the Sgr B2 (M) and (N) area, but this may be due to absorption, 
as shown in the CCS 45.38 GHz line (Fig©. This is different to 
th e distribution of the CC S 7,6 - 6,5 (81.505 GHz) line shown 
bv lKuan & Snvden 1 Jl994h . at higher resolution, which peaks near 
Sgr B2 (M) and (S). 



3.11 CH 3 CN 

There is one line of methyl cyanide CH3CN, shown in Fig. [5] the 
2(0) - 1(0) group (36.80 GHz). The emission also peaks in the north 
cloud, and shows the ridge in the integrated emission to the north- 



west of Sgr B2 (N) and (M), similar to that i n the CH3CN 5 - 4 
(91.99 GHz) and 6 - 5 (110.38 GHz) lines inljones et alj J200& 
and the 5-4 line in lde Vic ente. Martin-Pintado & Wilson d 199 
There may also be some absorption in the 2(0) - 1(0) line here, 
around Sgr B2 (N) and (M), not just a deficit of emission. 



3.12 OCS 

There are two lines of carbonyl sulphide OCS shown in Fig. [6] 
the 3 - 2 (36.49 GHz) and 4 - 3 (48.65 GHz) lines. The emission 
peaks in the north cloud, with mean peak position I — 0.684 deg, 
6 = —0.028 deg, mean velocity 63 km s _1 and mean velocity 
width 22 km s — , although there are some differences in distribu- 
tion between the two transitions. The distribution is similar to the 
7-6 (85.14 GHz) 8-7 (97.30 GHz) and 9 - 8 (109.46 GHz) lines 
(Jones et al. 2008a), with the ridge near Sgr B2 (N) and (M). 
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Figure 5. The integrated emission images for one line each of SO, SiO, 2J SiO, CS, ld CS and C a4 S, two lines of CCS and one line of CH3CN. Note that the 
dashed contours indicate negative integrated emission, or absorption. 



3.13 HNCO and HOCO+ 

There is one line of isocyanic acid HNCO the 2(0,2) - 1(0,1) group 
(43.96 GHz) and one of protonated C0 2 HOCO+ 2(0,2) - 1(0,1) 
(42.77 GHz) shown in Fig. [6] These lines peak at the north cloud, 
so the two molecules are discussed here together. 

The HNCO 43.96 GHz distribution is similar to that of the 
4( 0,4) - 3(0,3) (87.93 GHz ) and 5(0,5) -4(0, 4) (109.91 GHz) lines 
in I Jones et all d2008al) and lMinh et alj Jl998h . Observations of the 
4(0,4) - 3(0, 3) line over a wider a rea (albeit at lower resolution) 
are shown in Minh & Irvine] |2006) along with discussion of this 
north cloud (called the 2'N HNCO peak by them). T here are also 
observ ations of the 1(0,1) - 0(0,0) (21.98 GHz) line in lWilson et ail 
dl996l) . with similar distribution. 

The HOCO+ 42.77 GHz distribution similar to that of 
the 4(0,4) - 3(0 3) (85.53 G Hz) an d 5(0,5) - 4 (0 4) (1 06.91 
GHz) lines in Ijones et al.1 d2008al) . iMinh et alj d 19981) and 



Minh, Irvine & ZiurvsNl988l) and similar to the HNC O lines, high- 
lighti ng the chemical distinctness of this north cloud dwilson et al.1 
1996) from the Sgr B2 (M) and (N) area. 



3.14 CH 2 NH, NH 2 CN, HCS+ and CH 2 CN 

Finally, there is one line of each of the methylenimine CH2NH 
3(0,3) - 2(1,2) group (35.07 GHz), cyanamide NH 2 CN 2(1,2) - 
l(l,l),v=0 (39.73 GHz), thioformyl HCS+ 1 - (42.67 GHz) and 
the cyanomethyl radical CH 2 CN 2-1 group (40.25 GHz) shown 
in Fig. [6] 

The CH 2 NH 35.07 GHz line distribution peaks near 
Sgr B2 (M), and is simila r to the dis t ributio n of the 4(0,4) - 
3(1,3) (105.79 GHz) line in I Jones ^t aD |2008a). The latter being 
at higher resolution resolves the structure better, showing the peak 
at Sgr B2 (N) whereas the fitted peak in the 35.07 GHz line is near 
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Figure 6. The integrated emission images for two lines of OCS, and one line each of HNCO, HOCO+, CH 2 NH, NH 2 CN, HCS+ and CH 2 CN. 



Sgr B2 (M). There are likely to be excitation differences between 
the two cores. 

The NH 2 CN 39.73 GHz line peaks at the north cloud, with the 
ridge to the north-west of Sgr B2 (M), with a d istribution similar t o 
that of the 5(1,4) - 4(1,3) (100.63 GHz) line jjones et alj|2008ah . 
although the image there is noisy and affected by scanning stripes. 

The HCS + 42.67 GHz line is quite weak, but has an extended 
distribution. Not much is known on the distribution of this ion in 
Sgr B2, but it is plausible that the distribution is similar to that of 
CS (e.g. Fig.[5] subsection |3~9l. 

The distribution of integrated emission of the CH 2 CN 40.25 
GHz line is extended, with likely absorption near Sgr B2 (N). There 
are some stripes in the integrated emission in the longitude scan- 
ning direction. The distribution differs from that of CH3CN (Fig. 
[5] subsection |3.1U which may indicate that the two mo lecules have 
quite different formation routes and destruction paths l llrvine et"al] 
ll988l ; rTurneretal.lll990h despite their apparent similarity. 



3.15 Summary of peak positions 



The peak positions of the 7-mm lines in the data cubes tabulated 
in Table|2]are plotted in Fig. [7] The hydrogen recombination lines 
(filled triangles in Fig. [7} peak between Sgr B2 (N) and (M), as 
discussed in subsection 13.11 and shown in Fig. [2] The molecular 
lines mostly peak at the north cloud. 



We attribute the molecular line distributions in this 7-mm band 
to a combination of the effect of emission, from an extended cloud 
centred to the east of Sgr B2 (N) and (M), and absorption due to the 
free-free continuum at Sgr B2 (N) and (M). We consider this further 
in the Discussion (Section [5} and in the next section (Section [4]l 
where we analyse the spatial distributions in a relatively objective 
way. 
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Figure 7. The positions of the fitted 7-mm peaks, with filled circles for the 
molecular peaks and filled triangles for the H radio recombination lines. 
The contours are from the CH3CN line. The crosses are radio peaks and 
open squares infrared peaks, as in other images. 



4 PRINCIPAL COMPONENT ANALYSIS 

The integrated emission images of the different lines shown in Figs. 
[TJand[3]to[6]show similarities and differences. One useful technique 
to identify and quantify the features o f large data sets is Princi- 
pal Component Analysis (PCA), see e.g lHever & Schloerbl Jl997h . 
This describes the multi-dimensional data set by linear combina- 
tions of the data that describe the largest variance (the most signifi- 
cant common feature) and successively smaller variances (the next 
most significant features). 

In the context here, we can use PCA to describe the large num- 
ber of images, with a smaller set of image s which contain the mos t 
significant features. This has been used by Unge rechts et alj J 19971) 
for the OM C-1 ridge, an d more recently by us for the G333 molecu- 
lar cloud ( Lo et al. 2009) and the Sgr B2 area with 3-mm molecular 
lines jjones et alj|2008ch . 

We have implemented the PCA processing in a PYTHON 
script, with the PCA modulfl and pyFIT^I to read and write the 
FITS images. 

The images of the three most significant components of the 
7-mm integrated emission are shown in Fig. [8] The first four com- 
ponents describe respectively 56, 14, 9 and 4 percent of the vari- 
ance in the data. These components are statistical descriptions of 
the integrated line images, not necessarily physical components of 
Sgr B2. However, they do highlight the physical features in a useful 
way. 

The first component (Fig. [8} highlights emission from the 
north cloud and ridge. The second component highlights absorption 
(or emission) from Sgr B2 (N) and (M), and also some differences 
in the north cloud. The third component highlights differences in 



8 http://folk.uio.no/henninri/pca_module/ 

9 http://www.stsci.edu/resources/software_hardware/pyfits/ 



emission to the south. The fourth component (not shown here in 
Fig.[8j highlights differences in emission from Sgr B2 (N). 

The relations between the molecular integrated emission im- 
ages, and the PCA images are shown in Fig. [9] as vectors of the 
projection of the data images on the axes of the PCA images. The 
molecules are labelled with numbers, for clarity, in these plots. All 
of the data images are positively correlated with the first PCA im- 
age, except for SO (37). The greatest positive correlation with the 
second PCA component are for SO (37) and SiO (36), which have 
strong absorption at Sgr B2 (N) and (M). The greatest negative cor- 
relation with the second PCA component are for CH2NH (6) and 
the CH3OH 44.07 GHz line (12, with some maser emission), which 
have emission at Sgr B2 (N) and (M). The greatest positive corre- 
lation with the third PCA component are for the CH3CHO 38.512 
GHz line (7) and the H 13 CCCN 35.27 GHz line (17), and greatest 
negative correlation with the third PCA component are for NH2CN 
(33) and CS (15), highlighting more or less emission to the south, 
relative to the other lines. 

We used the integrated emission images for 37 molecules here 
in the PCA, excluding CH2CN (as it is probably affected by the 
scanning ripples), and not including the radio recombination lines 
(RRLs), as they are due to a different physical mechanism. How- 
ever, we have also calculated the PCA for the 46 integrated emis- 
sion images, including the RRLs, and find very similar results. The 
first four components then describe respectively 49, 24, 8 and 3 per- 
cent of the variance in the data, with the PCA images quantitatively 
very similar, with the major difference being the second component 
image having a stronger peak at Sgr B2 (N) and (M), describing the 
RRLs, and positive rather than negative. However, the sign of the 
component images is not physically significant. 

The PCA images from these Sgr B2 7-mm lines can be com- 
pared t o those fou nd for the Sgr B2 3-mm lines from I Jones et al.1 
.2008eJ) in Fig. [TO ( [jon"es et alj2008cb . The 3-mm PCA images are 
at higher resolution. 

The first two PCA component images are similar: the first 
highlighting the ridge and north cloud, and the second highlight- 
ing absorption at Sgr B2 (N) and (M). The second PCA component 
at 3-mm also includes emission to the south-west, as some of the 
strongest lines at 3-mm, such as HCN, HCO + and HNC are op- 
tically thick with both absorption at Sgr B2 (N) and (M) and the 
south-west area relativly stronger (as the the main emission is rela- 
tively weaker due to the high optical depth). The third 3-mm PCA 
component highlights the north cloud, and Sgr B2 (M). As PCA is 
a statistical tool, it is not expected that the higher PCA images, of 
lower statistical weight, should be the same from the independent 
7-mm and 3-mm data. However, the PCA images from the 3-mm 
and 7-mm data do show quite similar spatial features, albeit dis- 
tributed differently amongst the PCA images: ridge, north cloud, 
Sgr B2 (N) and Sgr B2 (M). 



5 DISCUSSION 

We can consider the distribution of the spectral lines in terms of 
physical components of the Sgr B2 complex: the hot cores Sgr B2 
(N), (M) and (S), and the extended envelope. These are related to 
the Principal Components discussed in the last section, but are not 
the sam e decomposition of structure into features. In Ijones et al.1 
(2008a) we identified seven features. Of these, the extended ridge 
and north cloud ar e related to the extended envelope of Sgr B2 
Jjones et al]|2008al) which is traced by the dust seen at sub-mm 
wavelengths (Fig. [A]}. The hot cores Sgr B2 (N), (M) and (S) are 
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Figure 9. The component vectors of the integrated line images in the decomposition of the data into the PCA images, shown with successive pairs of PCA 
images. The vectors are labelled with numbers, rather than molecule names, for clarity of presentation, with the molecules listed in alphabetical order: 1, 13 CS; 
2, 29 SiO;3, C 34 S; 4 and 5, two lines of CCS; 6, CH 2 NH; 7 to 10, four lines of CH 3 CHO; 11, CH3CN; 12 to 14, three lines of CH3OH; 15, CS; 16 and 17, 
two lines of H 13 CCCN; 18 to 24, seven lines of HC5N; 25, HCC 13 CN; 26 and 27, two lines of HC 3 N; 28, HCS+; 29, HNCO; 30, HOCO+; 31 and 32, two 
lines of NH 2 CHO; 33, NH 2 CN; 34 and 35, two lines of OCS; 36, SiO, and; 37, SO. 




GLON (degrees) GLON (degrees) GLON (degrees) 

Figure 10. The fir st three images from the Principal Component Analysis (PCA) for the 3-mm lines (lower) from ljones et al l J200 8ah. These 3-mm images are 
similar to those in ljones et alj j2008ct) rotated into Galactic coordinates for easier comparison to the 7-mm images. The first image highlights emission from 
the north cloud and ridge. The second image highlights absorption (or emission) from Sgr B2 (N) and (M) and emission from the south-west. The third image 
highlights the north cloud and differences in emission from Sgr B2 (M). 
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denser and warmer, and so have different chemistry, so they are 
prominent in particular molecules, as shown in Fig.[6]by CH2NH. 
In the 7-mm band data here, there is also often absorption in the 
area around Sgr B2 (N) and (M) related to the radiative transfer of 
the continuum free-free emission through the molecular gas. 

To better describe the quantitative properties of these physi- 
cal components of Sgr B2, we have considered in the Appendix 
continuum data (imaging and photometry) from the radio, through 
sub-mm to mid-IR wavelenths, from the literature, plus some of 
our own (as yet unpublished) radio and sub-mm data, and data 
from public archives. We have fitted the spectral energy distribu- 
tion (SED) of the Sgr B2 (N), (M) and (S) cores, and the extended 
envelope, using the imaging to set the angular size of the compo- 
nents, to derive parameters such as dust temperature. 

The integrated line emission J Tsdv plotted in Figs.[3]to[6]is 
related the the column density of the molecule in the upper state of 
the transition N u by 

N u = {8nv 2 k/hc 3 A ul ) J T B dv (1) 

where A u i is the Einstein coefficient, and assuming optically thin 
emission (optical depth r <C 1). This is related to the total column 
density of the molecule N, assuming local thermodynamic equilib- 
rium (LTE) by 

N = N u (Q T /gu) exp(E u /kT ex ) (2) 

where g u is the statistical weight of the upper level, Qt is the par- 
tition function at excitation temperature T ex and E u is the energy 
of the upper level. This analysis also assumes the line emission fills 
the beam (filling factor unity) and that the line brightness tempera- 
ture Tb is much greater than the continuum brightness temperature 
Ta- 
in this simple case, then, the integrated line emission J Tsdv 
traces the molecule column density N, and we can use this to trace 
chemical variations in the Sgr B2 complex. However, the simplify- 
ing assumptions are often not valid, so the relation is more compli- 
cated. 

Firstly, some transitions are not in LTE, as shown by the differ- 
ences between the three methanol (CH3OH) lines in Fig. [4] which 
include maser transitions. Secondly, even if the lines are in local 
thermodynamic equilibrium, there are large gradients in the kinetic 
temperature Tk and radiation environment (corresponding to ra- 
diation temperature T r ) giving gradients in excitation temperature 
T ex . So transitions from higher energy states will preferentially 
highlight higher temperature regions. Thirdly, the lines are often 
not optically thin, so that the optically thick emission is not simply 
proportional to the column density N u , in the sense that the inte- 
grated intensity saturates for high optical depth in regions of high 
column density. Fourthly, we have, so far, ignored the full radia- 
tive transfer model in which background continuum radiation can 
be absorbed by the molecules in the cloud, as well as the molecules 
emitting. With a background continuum source of brightness tem- 
perature Tc, we obtain 

T B = T ex [l - exp(-r)] + T c exp(-r) 

= T ex + (T c - T ex ) exp(-r) (3) 
« T c + (Tex-Tc)r for r < 1. 

Typically, the spectral line data (as in this paper) are baselevel sub- 
tracted, to be the line only, but the effect of the background con- 
tinuum is to reduce the line intensity, and for Tc > T ex the line 
becomes negative in intensity, i.e. is seen in absorption. 

Most of the lines presented here are dominated by emission 



from the extended, cooler (Tdust ~ 20 K), envelope, rather than the 
compact hot cores. The scale of this envelope is around 2 arcmin, 
but elongated north-south in equatorial coordinates (as shown in 
the mid-IR, Fig. lAU . This envelope is centred to the east of the line 
of hot cores Sgr B2 (N), (M) and (S), with the north cloud to the 
north of the hot cores. In Galactic coordinates, this corresponds to 
a ridge, wrapping around the hot cores, as shown in the first prin- 
cipal components of the 7-mm (Fig. [8} and 3-mm lines (Fig. 1101 ). 
There are also chemical differences, for example with HNCO and 
HOCO + preferentially found (Fig.[6]l in the north cloud. The effec- 
tive spatial scale of the stronger, optically thick lines, such as SO, 
SiO and CS (Fig. [5} is larger, with the emission filling the whole 
6x6 arcmin 2 area imaged, as the central peak emission saturates 
(and there is absorption). The optically thin isotopologues ( 29 SiO, 
13 CSandC 34 S, Fig.O are better tracers of the extent of the enve- 
lope. 

The recombination lines (Fig. |T} show the distribution of the 
7-mm free-free radio continuum, assuming that there is not a very 
large spatial variation in the line to continuum ratio (so that the 
line emission traces the continuum distribution). The deconvolved 
angular size (from these recombination lines) is 2.2 x 0.9 arcmin 2 , 
and the fit to the radio flux densities (Fig. I A2 1 and Table lA2t shows 
that around half the flux comes from the extended envelope and half 
from the compact hot cores. The spatial distribution along the line 
of sight of the Sgr B2 components is not well constrained, but the 
spectral energy distribution (SED) fitting (Appendix) does indicate 
that Sgr B2 (N) is behind at least part of the dust envelope. Many of 
the 7-mm molecular lines in Figs. [3] to [6] show absorption centred 
near Sgr B2 (N) and (M) which matches well the position and scale 
of the 7-mm free-free continuum (as probed by the recombination 
lines in Fig.[TJl. 

With a flat spectral index (a = -0.17± 0.07 for the total free- 
free) in the optically thin regime, the continuum brightness tem- 
perature falls steeply with frequency Tc oc v a ~ 2 . The conditions 
for absorption of the continuum radiation, by the intervening cool 
envelope (Tc significant compared to T ex ), are commonly met in 
these 7-mm data, as highlighted in the second principal component 
of Fig. [8] The molecules with the most significant absorption are 
SO, SiO, 29 SiO, CS, CH3CHO, NH2CHO and CCS (although for 
the last three molecules this is less clear in some of the lower signal 

to noise lines). 

Similar absorption occurs for the 3-mm data of Jones et al. 

(2008a) (seen in the second principal component of Fig. I10t . At 
these higher frequencies, the absorption at the hot cores Sgr B2 (N) 
and (M) dominates over the absorption due to the extended enve- 
lope. The resolution is higher, the cores contribute a larger fraction 
of the free-free continuum, and have a higher brightness tempera- 
ture. 

In both the 7-mm and 3-mm molecular lines, there are also 
lines that peak at the Sgr B2 (N) and (M) hot cores, due to their 
different chemical conditions and higher temperature. In the 7-mm 
data here, the most prominent 'hot core' lines are from CH2NH, 
and the 4 4.07 GHz transition of CH 3OH which includes maser 
emission (Mehringer & Mentenll 19970 . 

The chemistry of molecular clouds is complex 

( Herbs t & van Dishoecld [2009), particularly those associated 
with star-formation, like Sgr B2. Broadly, the chemistry of regions 
that are heated by star formation to become hot cores can be 
described in three stages: the cold, warm up and hot core stages. 

In the cold stage, there are gas-phase reactions plus grain- 
surface reactions, with a lot of the complex chemistry occuring be- 
tween molecules sticking on the ice-grains. In the warm up stage, as 
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the material forming the protostar collapes, there is heating (from 
10 K to a few hundred K), and the grain mantles sublime, releas- 
ing different molecules at different temperatures and hence times. 
There are further gas-phase and grain-surface reactions occuring 
with modified chemical abundances and hotter temperatures. In the 
last hot core stage, with temperatures 100 to 300 K, the ices have 
sublimed, so there is only gas-phase chemistry. 

The overall chemical evolution is complicated, and time de- 
pendent, so detailed models typically consider a particular phase. 
A recent gas-grain warm up mod el, of relevance to Sgr B2 (N) is 
iGarrod, Weaver. & Herbstl fl2008l) . 

The observed spectral-line data can be compared to such 
models, and constrain model s, but this is not straight-forward 
l lHerbst & van Dis hoeck 200jJ. For a particular position, or an un- 
resolved source, the observed antenna temperature can be con- 
verted to a source brightness temperature, using the source size or 
filling factor. This can be converted to a molecule column density 
(or abundance relative to e.g. CO) using the excitation temperature 
T ex and optical depth r (preferably using constraints from other 
lines). A better way to analyse the molec ular line data is to con- 
struct a physical model of the source (see iHerbst & van Dishoeck! 
120091 Figure 4) and calculate the full non-LTE statistical equilib- 
rium excitation of the molecule through the physical structure. 

For the Sgr B2 complex data presented here, with the hot 
cores, and cooler envelope, we do not attempt such a detailed radia- 
tive transfer model, but simply point out the broad features required 
in such a model, namely the physical components (see Appendix) 
and the effects of optical depth and absorption. We show by the 
line distribution images that there are significant spatial differences 
in the chemical distribution, that can be interpreted in terms of the 
cool, warm up and hot core 'stages' of the chemical models, with 
different parts of the Sgr B2 complex having different heating con- 
ditions and history. 



6 SUMMARY 

We have undertaken a line imaging survey of the Sgr B2 region 
from 30 to 50 GHz (the 7-mm band) with the 22-m single-dish 
Mopra telescope. This co mplements t he Mo pra 3-mm line imag- 
ing of Sgr B2 presented in lJones et alj fl2008ah . Integrated emission 
images of 47 lines are presented: 38 molecular lines and 9 radio 
recombination lines. 

The distributions of individual lines are discussed, and we 
have studied the similarities and differences between the lines with 
principal component analysis (PCA). The major features from the 
PCA are extended emission from the ridge and north cloud to the 
north and east of the Sgr B2 (N) and (M) hot cores, and absorption 
near Sgr B2 (N) and (M). These statistical features are interpreted 
in terms of the physical components of the Sgr B2 complex (dis- 
cussed in the appendix) of the extended low temperature envelope, 
and compact hot cores Sgr B2 (N), (M) and (S). The 7-mm free- 
free emission in the area of Sgr B2 (N) and (M) is absorbed in this 
line-of-sight by the envelope. 
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APPENDIX A: THE COMPONENTS OF SGR B2 

The Sgr B2 complex is very complicated, with many physical 
components distingui s hable at hi gh resolution , see f or example 
iMehringer & Mentenl d 19971) and fcaiime et al] dl995h . However, 
for many observations with resolution of order tens of arcseconds 
(Fig. lAlt it is common, and physically useful, to consider the three 
dense cores Sgr B2 (N), (M) and (S), and distinguish them from the 
larger and less dense envelope. 

The three dense cores are prominent in both the radio free- 
free and sub-mm dust continuum, as shown in Fig. lAll and the (M) 
and (S) cores are prominent in the mid-infrared. The sub-mm dust 
distribution also shows the extended component, centred slightly to 
the west of the line of dense cores, wit h an extension to the north. 
As we pointed out in ljones et al. (2008a), many of the 3-mm molec- 
ular lines follow the ridge of this extended dust component, while 
other lines are associated with the dense cores. We find similar dis- 
tributions here for the 7-mm molecular lines, so for more detailed 
quantitative discussion of these 7-mm and 3-mm line data, we have 
undertaken an analysis of the Sgr B2 continuum data. 

The structure o f the Sgr B2 complex has been analysed by 
Gord on et "all dl993h using the spectral energy distribution (SED) 
of the total complex and components FIR1, FIR2 and FIR3 in their 
nomenclature, corresponding to (N), (M) and (S). 

The peak in the SED due to dust is modelled as a 
Black Body modified by the optical depth r„ of the dust, 
so the flux density S v divided by the source solid angle Q. 
is S v /n = B v (T)(l - exp(-T„)) where B V (T) is the 
Planck function (Wilson, Rohlfs, Hiittemeister 2009). The dust op- 
tical depth is assumed to be proportional to the gas column 
density, with some empirica l scaling factor. The model from 
iMezger. Zylka & Winkl jl990l) is used with r„ = (/Vsr/1.4 x 
10 20 cm" 2 )(A/Atm)" 2 6(Z/Z Q ), where N H is the total hydro- 
gen column density, (Z/Zq) is the metallicity, assumed 2 for the 
Galactic Centre and b is a dust parameter, with b = 3.4 for hydro- 
gen number density nn 10 6 cm -2 and b — 1.9 for nu ^ 10 6 
cm~ 2 . 

The observed flux densities in the dust SED can be fitted with 
3 parameters; temperature T, solid angle and the wavelength 
A r= i at which the optical depth is unity, related by the above equa- 
tion to the column density Nh ■ (In this model the power-law index 
of the dust optical depth is /3 = 2: other models can have different 
power law slopes /?, but we leave /3 fixed here to reduce the number 
of free parameters in the fit.) 

The SED model assumed is simple, with constant temperature 
T, over a uniform dust density distribution of solid angle fi, but 
this is a common assumption in fitting sub-mm dust SEDs. It also 
has the advantage of a simple analytic form that can be fitted. We 
recognise that the dust distribution and temperature are not uniform 
(as shown in higher resolution studies) but argue that the character- 
istic temperature, density and other parameters derived (Table [A2t 
are useful for quantitative analysis. 

We plot in Fig. I A2 1 the SED of the Sgr B2 components, using 
the data from lGordon et al] d 19931) plus some additional data, from 
more recent observations. 

We include fluxes from lKuan. Mehrin ger & Snvdeil d 19961) for 
the compact components in the 3-cm and 3-mm bands (combining 
their N and N' into the one source). We also include mid-infrared 
fluxes from the MSX catalogue Baby Gator interfacj"! 

We have fitted data from our own observations with the ATCA 



http://irsa.ipac.caltech.edu/applications/BabyGator/ 
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RA (J3000) RA (J3000) RA (J2000) 

Figure Al. The continuum images of the Sgr B2 complex in radio (ATCA 1.54 GHz), sub-mm (SCUBA 850 (im) and mid-infrared (MSX 21 £im), showing 
the Sgr B2 (N), (M) and (N) cores. The radio and mid-infrared trace star formation, with some extra compact features, and extended emission to the south-east. 
The sub-mm traces warm and cold dust, including diffuse emission around the cores. Note that these images are in equatorial coordinates, rather than Galactic 
coordinates used in the main body of this paper, but the crosses for radio peaks and open squares for infrared peaks are used as fiducial marks. 



log(Wavelength / micron) 

_l I I I L_ 



log(Wavelength / micron) 




12 3 4 
log(Frequency / GHz) 

log(Wavelength / micron) 





12 3 4 
log(Frequency / GHz) 

log(Wavelength / micron) 




12 3 4 
log(Frequency / GHz) 



1 2 3 4 5 
log(Frequency / GHz) 



Figure A2. The spectral energy distribution (SED) of the Sgr B2 components (N), (M), (S) and total, from the radio to the infrared, following Gord on et alj 
1 1993), with some additional data. The dashed lines are our schematic fits to the SED, but we note that there are problems using fluxes from the literature with 
a wide range of spatial resolutions so that the x 2 value of the model fits are poor. The plot for the total flux has a long dashed line for the total flux points, and 
a shorter dashed line for the extended envelope which contains most of the flux. The plot for Sgr B2 (N) includes 3-sigma upper limits, from Spitzer MIPS at 
70 and 24 micron, and MSX at 21 micron. 
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Table Al. Fitted flux densities used for the Sgr B2 spectral energy distribu- 
tion, in addition to values from the literature. 



Frequency 


Telescope 




Flux density 




or 


Project or 


Total 


N 


M 


S 


Wavelength 


Instrument 


Jy 


Jy 


Jy 


Jy 




ATCA 




2.35 


2.57 


13 


y. 1 Ljrlz 


VLA 




6.3 


7.9 


U.o4 


JO Vrn/ 


ATC A 




6.1 


8.0 




1 .2 mm 


MMdA 


tut 


103 


7 1 


J. 1 


1.1 mm 


BOLOCAM 


420 


64 


90 


25 


0.85 mm 


SCUBA 


2500 


330 


260 


22 


0.45 mm 


SCUBA 


68000 


3200 


3300 


720 


100 itm 


IRAS 


65000 








60 fim 


IRAS 


14500 








70 fim 


Spitzer 




<280 






24 ttm 


Spitzer 




< 0.40 






21 /im 


MSX 




<2.1 


16.4 


9.5 


15 fj,m 


MSX 






1.80 


2.13 


12 Atm 


MSX 






1.07 


0.96 


8 Aim 


MSX 






0.79 


0.83 



at 1.54 GHz (Fig. [AlJ and 36 GHz (Fig. [2jl and with SIMBA on 
SEST at 1.2 mm. The fitting was done with multiple gaussians us- 
ing the task imfit in the MIRIAD package, to obtain angular size 
information as well as flux densities. 

We have also fitted data obtained as FITS images from sev- 
eral public databases. In the radio we obtained VLA images at 20 
cm, 3.6 cm, 9 GHz (3.3 cm) and BIMA data at 3 mm from the 
Astronomy Digital Image Library (ADIL0 In the sub-mm we ob- 
tained JCMT SCUB^3 images at 450 and 850 Atm an d recent CSO 
BOLOCA1V0 images at 1.1 mm dBallv et aljHoToh . In the mid- 
infrared we obtained MS^f^l images at 8, 12, 15 and 21 fj,m, Spitzer 
MIPS[3 images at 24 and 70 /jm and IRAS reprocessed HIRES'^ 
images at 60 and 100 /im. The recent paper (Bally et alj|2010h dis- 
cusses the BOLOCAM 1.1 mm data, and shows the Spitzer 24 Aim 
as a greyscale in their Figure 18 with 1.1 mm contours. They also 
show 350 fim CSO SHARC-II data, but do not quote fluxes (and 
we did not have access to the 350 Aim from public archives). 

The extra flux density data are summarised in Table I All The 
limits for Sgr B2 (N) at 70, 24 and 21 Aim are 3-sigma upper limits 
(see later discussion). 

The flux density data (Fig. IA2I > show two physical compo- 
nents, the dust peak around 1.5 THz (200 fim) and the optically 
thin free-free emission between 3 and 100 GHz. There are, how- 
ever, problems in using flux densities obtained from data obtained 
with a range of spatial resolutions. For the fitting of the complex 
structure with the components Sgr B2 (N), (M) and (N), and to- 
tal flux, different amounts of flux will be included depending on 
the size of the telescope beam (as the four components are a sim- 
plification). This leads to discrepancies in the flux densities, even 



11 http://imagelib.ncsa.uiuc.edu/imagelib.html. Note that some of the 
archived radio images from ADIL were not suitable for measuring flux den- 
sities, as the images were not primary beam corrected, but they were useful 
for measuring sizes of the compact components. 

12 http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/jcmt/ 

13 http://ksa.ipac.caltech.edu/data/BOLOCAM.GPS/ 

14 http://irsa.ipac.caltech.edu/applications/MSX/MSX/ 

15 http://irsa.ipac.caltech.edu/applications/Spitzer/Spitzer/ 

16 http://irsa.ipac.caltech.edu/IRASdocs/hires_over.html 



considering the quoted uncertainties. We have fitted models to the 
SEDs, despite this, but note that the reduced \ 2 of the fits are high. 

For the free-free component, we have fitted power laws over 
the frequency range up to around 100 GHz, excluding some points 
for Sgr B2 (N) and the total which show some upturn due to the 
dust component. For the dust component we fitted the 3-parameter 
model disc ussed above. For Sgr B2 (N), as we discuss below, we 
also follow iGordon et al.1 d 19931) in including an extra absorption 
term for dust along the line of sight, making a 4-parameter fit 
S v = Q.B V (T)(1 — exp(— t„.i)) exp(— 7V,2). The dust fit used 
the data around 100 GHz for Sgr B2 (N) and the total, after cor- 
rection for the free-free component. The fitting was done in the 
GNUPLOT package, which uses the nonlinear least-squares (NLLS) 
Marquardt-Levenberg algorithm. For the total flux, we fitted the 
data after subtracting the models for the compact components Sgr 
B2 (N), (M) and (S), so that we were actually fitting the model to 
the extended envelope. The model fits for both the total flux and the 
extended envelope are both plotted in Fig. lA2l and we note that for 
most of the frequency range, the extended component dominates 
the total flux. 

For the free-free component, we find spectral indices of 
-0.02 ± 0.03, 0.22 ± 0.02, 0.04 ± 0.16 and -0.32 ± 0.12 for 
the components Sgr B2 (N), (M), (S) and the extended envelope 
respectively, and flux densities at 30 GHz 7.1 Jy, 9.9 Jy, 0.62 Jy 
and 17.7 Jy respectively. Note that this indicates that around half 
the flux is in the extended component which may be missed by 
the interferometric observations (e.g. Fig [2jl. We do attribute the 
radio emission from 3 to 100 GHz in the compact components to 
optically thin free-free emission, with basically a flat spectrum, al- 
though the rising spectrum of Sgr B2 (M) may indicate some of 
the subcomponents are ultracompact H II regions with the turnover 
from optically thick to optically thin above 10 GHz. The spectrum 
of the total Sgr B2 flux is flat (spectral index -0.17 ± 0.07), but 
when we subtract the compact components, notably Sgr B2 (M) 
with the rising spectrum, the residual extended component has a 
somewhat steeper spectrum. 

There has been s ome discussion of n o n-thermal radi o emis- 
sion from Sgr B2 by iHollis et al.1 ( 120071) . ICrocker et al.1 d2007l) 
and IProtheroe et al.1 d2008h. We agree with lLang. Palmer & Gossl 
(2008) and lProtheroe et al.l l2008) that due to the complex structure 
of Sgr B2, it is difficult to measure the radio spectral index accu- 
rately by combining observations with different resolutions which 
may be including diffe rent features . In pa rticular, the fluxes mea- 
sured with the GBT bv lHollis et all d2007t) for Sgr B2 (N) will in- 
clude flux from the extended component, to a greater extent at the 
lower frequencies with the larger beam. 

The fits to the dust components, for thr ee free param e ters T , 
Q and A r= i gave similar results to that of IGordon et all d 1993b . 
as expected since the data used here is based on their compilation. 
However, there are quite large uncertainties in derived Q and A T= i, 
and the two parameters are highly correlated. We decided, there- 
fore, to constrain the fits by using a fixed derived value for the dust 
solid angle Q obtained using the extra spatial information from the 
sub-mm images. This fit to two free parameters (Q and A T= i) gives 
similar reduced \ 2 to the three free parameter fits, and only changes 
the fitted temperature T by a few K. 

We used the geometric mean of the deconvolved sizes from the 
SCUBA 450- and 850-Aim images (resolution 7.5 and 14 arcsec), 
and the 1.3-mm IRAM 30-m image (resolution 11 arcsec), giving 
gaussian half-widths 10.2, 11.2 and 10.2 arcsec for Sgr B2 (N), (M) 
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Table A2. Parameters for the Sgr B2 components fitted from the SEDs, and derived from these fits. From the fit to the optically-thin free-free part of each 
spectrum, we obtain the spectral index ctf f and the flux density S30 at reference frequency 30 GHz. We choose to fix the angular size of the dust component 
to that measured from the sub-millimetre images, and derive the dust temperature T d and turnover wavelength A T= i from the dust peak in the spectrum. The 
total hydrogen column density Nh, number density njj and mass .'./// are derived from the dust spectrum fits, as described in the text. 



Component 


a ff 


S30 


U 


Diameter 


r d 


A T= i 


N H 


n H 


M H 






fixed 








10 24 cm~ 2 


10 6 cm" 3 








Jy 


arcsec 2 


pc 


K 


/nn 


(10 28 m" 2 ) 


(10 12 m" 3 ) 


M 


N 


-0.02 ± 0.03 


7.1 


117 


0.40 


65 ± 17 


740 ± 150 


12 


8.9 


8600 


M 


0.22 ± 0.02 


9.9 


138 


0.43 


50 ± 1 


680 ± 70 


9.6 


6.8 


8400 


S 


0.04 ±0.16 


0.62 


118 


0.40 


48 ± 1 


230 ± 30 


1.1 


0.87 


850 


envelope 


-0.32 ±0.12 


17.7 


17400 


4.5 


24 ±2 


280 ± 40 


3.0 


0.19 


320000 



and (S).0We therefore used solid angles Q = 1.1336* 2 of 117, 
138 and 1 18 arcsec 2 for Sgr B2 (N), (M) and (S). For the extended 
component, we used the geometric mean of effective fitted solid 
angles from the SCUBA 450- and 850-^m, BOLOCAM 1.1 -mm 
and 1.3-mm IRAM-30-m images, which is Q. = 17400 arcsec 2 . 

The results from the dust spectrum fits, with fixed, and two 
free parameters T, and A T =i are given in Table lA2l 

For Sgr B2 (N) we find, as pointed out bv lGordon et al.l Jl993l) . 
that the extra line-of-sight absorption component to the model is 
needed, making in this case then a three parameter fit. If we fix 
the angular size at that measured, then the extra absorption com- 
ponent to the model is clearly required: the reduced x 2 of the fit 
without the absorption (7.7) is substantially greater than with the 
absorption included (4.4), and the fit is clearly poor for the sub-mm 
range. 

The turnover wavelen gth is fitted as A T =i = 100 ± 17/^m, in 
good agreement with lGordon et alj a 1993T ). This is physically inter- 
preted as Sgr B2 (N) being partly embedded in the envelope. The 
fit shown in Fig. IA2l does not use the upper limits as constraints, 
but the non-detection of Sgr B2 (N) in the Spitzer MIPS 24 micron 
data (at < 0.40 Jy) is consistent with the fit, and is a strong con- 
str aint. The fit does rely st rongly on t he 100 and 50 micron points 
of lGoldsmithetalJ dl992l) , quoted in iGordon etal] jl993h where 
Sgr B2 (N) was detected as a feature somewhat merged with the 
stronger peak of Sgr B2 (M) and the extended envelope. Our fitting 
of the recent Spitzer MIPS 70 micron image shows Sgr B2 (M) 
and (S) saturated, so that these components and the envelope can- 
not easily be fitted. The upper limit we quote for Sgr B2 (N) at 70 
micron is the flux at the known (from sub-mm) position, in addi- 
tion to so me extended structure, w hich indicates that the 50 micron 
flux from Golds mith et al.l Jl992h may be an over-estimate, includ- 
ing some flux from extended structure. This would indicate more 
absorption (and longer A T =i) may be appropriate, but better high 
resolution data in the far-infrared range are required to confirm this. 

Due to the inter-relationship of parameters T, £1 and A T= i it is 
possible to fit the cutoff in the dust SED at high frequencies (infra- 
red, where Sgr B2 (N) is not detected, see Fig. lAlt with angular size 
SI as a free parameter, which gives a low temperature (T around 30 
K), but this requires a large angular size S7 ~ 500 arcsec 2 which 
is in disagreement with the measured angular size. In addition, if 
Sgr B2 (N) had this low temperature, similar to that of the extended 



envelope, (and there was not the absorption), then the 70 micron 
Spitzer MIPS image would show a clear peak at the Sgr B2 (N) 
position, in addition to the extended structure. This is not the case, 
which is further evidence against the low temperature interpretation 
of the Sgr B2 (N) SED. 

We also list in Table lA2l the parameters derived from the dust 
spectrum fits of total hydrogen column density Nh, number den- 
sity uh and mass Mh for the four components Sgr B2 (N), (M), 
(S) and envelope. Note that N H = N(H I) + 2N(H 2 ) so this 
includes both molecular and atomic hydrogen. The turnover wave- 
length A T= i is related to total hydrogen column density Nh by 
Nh = 1.4 x lO 2O cm- 2 (A T =i/Vm) 2 (Z /Z)/b where Z/Zq 
is the metallicity, assumed 2 for the Galactic Centre, and b is 
a dust parameter, assumed 3.4 for components Sgr B2 (N), (M) 
and (S), with nn ^ 10~ 6 cm~ 6 and 1.9 for the envelope with 
n H s= 10~ 6 cm- 6 dGordon et al]|l993l ; iMezger et al.lll990h . The 
number density nn and mass Mh are derived using the angular 
size f2 to estimate the line-of-sight thickness, and hence volume of 
th e components. These are rough estimates, similar to those derived 
bv lGordon eFaT] ( 1 1993h from similar fits, but we assume Galactic 
Centre distance 8.0 kpc. 

We note that the three hot cores Sgr B2 (N), (M) and (S) have 
similar angular size and temperature, but that Sgr B2 (S) has around 
an order of magnitude smaller flux density (both dust and free-free 
components), column density, number density and mass. The ex- 
tended envelope is cooler but around an order of magnitude larger 
in diameter than the cores, or two orders of magnitude in solid an- 
gle, and dominates the flux density of both dust and free-free com- 
ponents. The envelope is of much lower number density than the 
cores, but with much greater mass. 

We use these parameters of the Sgr B2 components in Section 
[5]to interpret the molecular line observations of Section[3] 

This paper has been typeset from a TpX/ KTpX file prepared by the 
author. 



We did not include deconvolved sizes from the BOLOCAM 1.1-mm or 
SIMBA 1.2-mm images, as the images were lower resolution (33 and 24 
arcsec) and hence the deconvolved sizes were much less accurate. 
18 The SIMBA 1.2-mm image may suffer from spatial filtering in the data 
reduction, affecting the extended structure. 



© 2010 RAS, MNRAS 000.[Q[[8| 



